Excitation signal optimization for minimizing fluctuations in knock out slow extraction

The synchrotron is a circular particle accelerator used for high energy physics experiments, material and life science, as well as hadron cancer therapy. After acceleration to the desired energies, particle beams are commonly extracted from the synchrotron using the method of resonant slow extraction. The goal is to deliver a steady particle flux—referred to as spill—to experiments and treatment facilities over the course of seconds while slowly emptying the storage ring. Any uncontrolled intensity fluctuations in the spill are detrimental to the efficiency of beam usage, as they lead to detector pileups or detector interlocks, hindering experiments and cancer treatment. Among the most widely used extraction scheme in medical facilities is the Radio Frequency Knock Out (RF-KO) driven resonant slow extraction, where the stored beam is transversely excited with a radio frequency (RF) field and the spill intensity is controlled by the excitation signal strength. This article presents particle dynamics simulations of the RF-KO system with the focus on finding effective mechanism for minimizing the intensity fluctuations while maintaining a good extraction efficiency and other advantages of KO extraction. An improved beam excitation signal which optimizes these main objectives is found, and is rigorously compared experimentally with other commonly applied techniques.


Excitation signal optimization for minimizing fluctuations in knock out slow extraction Philipp Niedermayer * & Rahul Singh
The synchrotron is a circular particle accelerator used for high energy physics experiments, material and life science, as well as hadron cancer therapy.After acceleration to the desired energies, particle beams are commonly extracted from the synchrotron using the method of resonant slow extraction.The goal is to deliver a steady particle flux-referred to as spill-to experiments and treatment facilities over the course of seconds while slowly emptying the storage ring.Any uncontrolled intensity fluctuations in the spill are detrimental to the efficiency of beam usage, as they lead to detector pileups or detector interlocks, hindering experiments and cancer treatment.Among the most widely used extraction scheme in medical facilities is the Radio Frequency Knock Out (RF-KO) driven resonant slow extraction, where the stored beam is transversely excited with a radio frequency (RF) field and the spill intensity is controlled by the excitation signal strength.This article presents particle dynamics simulations of the RF-KO system with the focus on finding effective mechanism for minimizing the intensity fluctuations while maintaining a good extraction efficiency and other advantages of KO extraction.An improved beam excitation signal which optimizes these main objectives is found, and is rigorously compared experimentally with other commonly applied techniques.
Synchrotrons have proven to be vital instruments for studying the fundamental properties of nature.In the last two decades, they have also facilitated large scale effective hadron cancer therapy [1][2][3][4] .For slow beam extraction, the RF-KO method is implemented in most synchrotrons since it avoids any drift of the beam during extraction, enables a fast abort of the extraction and allows for spill intensity control [5][6][7] .From the user's perspective, the spill quality along with beam extraction efficiency are the most important figures of merit of a synchrotron facility.As any resonant slow extraction scheme, RF-KO extraction suffers from spill fluctuations caused by minuscule magnet power supply ripples, but in addition, artefacts of the RF excitation signals driving the extraction are imprinted on the spill.This is shown to cause severe impairments in effective beam usage, extensively reported in literature [8][9][10][11] .Numerous efforts have been made in the past seven decades to address this particular challenge [12][13][14][15][16][17][18][19] .
More recently, a multitude of heuristic improvements to the RF-KO process have been made in terms of tailored excitation signals [20][21][22][23] .This contribution reviews the RF-KO beam dynamics and presents a novel spill smoothing mechanism which is robust to the specificity of synchrotron settings and largely independent of any specific accelerator facility.The presented excitation signal is compared with all the other contemporary excitation signals by means of simulations and experiments.
After introducing the mechanism of RF-KO extraction along with the relevant beam dynamics and spill quality metrics in the rest of this section, the subsequent section focusses on the design of appropriate excitation signals along with simulation studies.A robust mechanism of spill smoothing is explained, leading to the proposal of a novel excitation signal referred to as Noise++.The following section experimentally compares the excitation signals under identical machine conditions and discusses the observed effects on the spill.

RF knock out resonant slow extraction
The general principle of resonant slow extraction is depicted in Fig. 1.The working point (tune) of the circular accelerator in the plane of extraction (typically horizontal) is moved close to a 3rd order resonance and nonlinear sextupole fields are used to drive the resonance 24 .The three-turn particle dynamic under these conditions is described by the Kobayashi Hamiltonian 25 (1)

OPEN
GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany.* email: p.niedermayer@gsi.de where S is the normalized sextupole strength and (X, X ′ ) are the normalized phase space coordinates with where the action J x and angle x were used and the relation to the physical coordinate x and divergence x ′ is given by the beta function β x (s) which describes the optical properties of the accelerator at the location s.The (small) quantity d = Q x − Q res is the distance of the tune Q x to the resonance Q res with 3Q res ∈ N .In the follow- ing, uppercase letters Q will be used to denote absolute tune values, while lowercase letters q ≤ 0.5 refer to the fractional tune as distance to the nearest integer.As a consequence of the non-linear dynamics in the vicinity of the resonance, a transverse instability is created in phase space (Fig. 1, left).At the stability limit (separatrix), the Hamiltonian takes the value H sep = (4π d) 3 /S 2 .When a particle transitions beyond the separatrix, its motion becomes unbound and its action increases exponentially.
For RF-KO driven resonant slow extraction, the particle action is controlled with transverse excitation 26 .For this purpose, an electromagnetic RF field created inside a kicker element (exciter) 27 deflects the traversing particles on each turn, leading to a successive increase of the action.After crossing the stability limit, the rapid, exponential increase allows particles to jump over the thin blade of an electric septum.The septum field deflects the particles such that a subsequent magnetic septum can be used to guide them into the extraction beam line (Fig. 1, right) 7 .This spill of particles is then delivered to experiments, patients or subsequent accelerators.

Spill fluctuations and quality
The number of particles being extracted per unit time (spill rate) fluctuates.Typically, these fluctuations are far more severe than one would expect from a Poisson process governing the discrete statistics of an uniform extraction rate of particles.In RF-KO extraction, these spill fluctuations have two primary causes, both of which are enhanced if the particle density in the vicinity of the separatrix is high.
The first cause are current ripples with relative magnitudes in the order of 10 −5 caused by the power sup- plies feeding the accelerator magnets 28 .During the conversion to direct current (DC), harmonics of the mains frequency (50 Hz)-which cannot be filtered out completely-remain as small ripples.Together with noise and ripples up to a few kHz stemming from control circuits in the power supplies, these perturb the magnetic guiding field, altering the tune and sextupole strength.As a result, the stability limit H sep ∝ d 3 /S 2 fluctuates, which leads to significant modulations of the spill rate during the extraction process 7 .
The second cause of spill fluctuation lies in the excitation driving the extraction process.The interaction of the excitation signal with the anharmonic oscillations of particles as their motion becomes increasingly non-linear in the vicinity of the separatrix affects the extraction process by modulating the particle coordinates and can introduce fluctuation in the spill rate.When using inappropriate excitation signals, this contribution is typically dominant.The design of excitation signals is therefore not only important to prevent such additional fluctuations, but also opens up the opportunity to further reduce spill fluctuations by using tailored excitation signals.
To quantify the spill fluctuations and provide a quality metric, the spill is divided into time intervals of length t count and the number of extracted particles N is counted in each interval.The choice of the interval length is subject to the temporal resolution of the detector system being used.The fluctuation of particle counts is evaluated over a larger time span t evaluate by means of the coefficient of variation c v or the related spill duty factor F which are defined as 28 (2) Circles (filled for every 3rd turn) show the trace of a particle crossing the separatrix (red dashed), becoming unstable and jumping over the septum blade (black).Right: Aerial view of the beamline.The exciter controls the transition from the stable (thin line) to the unstable orbit (thick line), where the particle is deflected by the electrostatic (E) and magnetic (M) septa.

Simulation framework
The particle tracking simulations are based on Xsuite 29,30 and use the lattice model of the Cooler Synchrotron (COSY) with the same beam and optics parameters as in the experimental investigation (Table 1).For computational reasons, the simulated time is reduced to 2 s ( 2.4 × 10 6 turns) and only 10 6 and 10 4 particles are tracked for the spill quality and beam dynamics studies, respectively.Power supply ripples are modelled by introducing a modulation of the quadrupole strength at 150, 300 and 600 Hz as well as noise below 10 kHz, motivated by the frequencies typically observed in spill spectra.This results in a periodic tune change of q ripple ≈ 10 −5 .
The excitation signals are generated with the GNU Radio software 31 .Inline with the experiments, a feedback controller adjusts the signal level to maintain the targeted spill rate, which is given by the number of simulated particles divided by the spill duration.The following signal parameters are used in the simulation: The excitation frequency f ex = Q ex f rev ≈ (n ± q x )f rev lies near a betatron sideband ( q x = 0.346 ) at the harmonic n ∈ N of the revolution frequency f rev .The broad noise band is centered on the first sideband at f ex = 0.346f rev = 414 kHz with a bandwidth of f ex = 0.028f rev = 33 kHz covering the full range to the resonance q res = 1/3 as traditionally performed.For the sinusoidal excitation a fixed frequency of f ex = 0.656f rev = 784 kHz is chosen as a result of a rough parameter scan for the lowest signal level required to maintain the extraction.Using a pure sinusoidal signal is only possible because of the non-linear dynamics involved, yet the root mean square (RMS) signal level required is about a factor 10 larger than for the noise signal.This aspect is further discussed in "Proposed (3) www.nature.com/scientificreports/excitation signal" section.The simulations are executed on a high-performance computing (HPC) cluster using dedicated graphics processing units (GPUs) and require computation times of several hours per spill up to a week.

Spill smoothing by coherent excitation
Figure 2 shows two simulated spills, one extracted with the sinusoidal signal and the other with the band limited noise signal.The essential aspect controlled by the excitation signal spectra is the time a particle spends in the vicinity of the separatrix, where it is subject to fluctuations.For the broadband excitation signal, the extraction is a slow diffusive process and a significant interaction occurs with the non-linear particle motion near the separatrix.This translates onto the spill by means of the low frequency characteristics visible in Fig. 2. The majority of this low frequency noise is not caused by the quadrupole ripples, but rather is a consequence of the diffusive extraction in combination with the typical low-pass behaviour of the extraction process 19 .In contrast to this, for sinusoidal excitation, the particle motion is coherent, leading to a fast crossing of the separatrix and minimal interaction with it.Thus, the imprinting of low frequency noise is avoided, and the influence of separatrix fluctuation due to magnetic field ripples is weakened.Consequently, the spill fluctuations are strongly reduced for coherent sinusoidal excitation compared to excitation with a non-zero bandwidth, resulting in a more uniform spill rate with fewer spikes.A detailed description of this process is provided in the following subsection.
An intrinsic property of the RF-KO extraction is the cause of coherent beam oscillations at the respective excitation frequencies, leading to an imprint of these frequencies and its harmonics on the spill spectrum (Fig. 2).The resulting peaks in the MHz domain are strongest for the fully coherent sinusoidal excitation, and they can be moved to higher frequencies if required, by using excitation frequencies corresponding to betatron sidebands at higher harmonics.For the considered spill rate of 5 × 10 5 particles/s the excitation frequencies used in this study do not impact the spill quality as they are sufficiently high.For higher spill intensities the periodicity f −1 ex of the excitation signal can be observed as periodic structure on the spill (Fig. 3).When the low frequency fluctuations described in the previous paragraph cause a momentary increase of the extraction rate, this increases the probability of multiple particles being extracted in short succession and causes pile-up.Although the sinusoidal excitation does not prevent these structures from being imprinted on the spill, the spill smoothing in the low-frequency domain avoids the occurrence of such short-term increases in the rate and thus reduces the probability of pile-up.At the same time, discontinuities in the spill with very low rate are avoided, ensuring an efficient use of the beam.

Non-linear particle dynamics under excitation
The reason for the reduced spill fluctuations can be understood by analysing the non-linear beam dynamics of the excited particle motion.Starting from the equations of motion defined by the Kobayashi Hamiltonian (equation 1), an additional term is derived, taking into account a small time-dependent sinusoidal excitation with a dipolar field modelled as a kick with the normalized dipole strength K proportional to the signal level, the turn number n = tf rev , the excitation tune Q ex and phase φ ex .After the kick, the coordinates of a particle initially at (X, X ′ ) read This corresponds to the action and angle (5) Figure 3. Spill rate on sub-µ s timescale from particle tracking simulation with an increased target rate of Ṅ = 10 7 /s for comparison of noise and sinusoidal excitation.The inset plot shows the shaded region with 100 ns resolution.
Since the kick X ′ is small, one can approximate In three consecutive turns, with the tune close to the 3rd order resonance Q res , the accumulated change is This additional change due to the small sinusoidal excitation can be expressed as a change H of the value of the Kobayashi Hamiltonian (equation 1) every three turns, such that Hamilton's equations �� x,3 = ∂�H/∂J x and �J x,3 = ∂�H/∂� x are fulfilled.This yields where n = 0, 3, 6, . . . is the turn number.
The sinusoidal excitation periodically increases and decreases the value of H by H every three turns, causing the particle to spiral inward and outward in transverse phase space as depicted in Fig. 4. The frequency of this modulation of H is determined by Q ex and � x ≈ 2πQ x n and is approximately |Q ex − Q x |f rev ≪ Q x f rev .Since x is itself modulated due to detuning, the frequency of this spiralling motion varies, especially for H → H sep where ( 8)  the detuning becomes stronger.The effect is also proportional to √ J x , meaning that it is particularly strong near the corners of the separatrix, i.e. the location where the particle can become resonant and unstable.Here, the excitation has the effect of pushing the particle periodically away from, or closer to the resonant case, eventually leading to a faster crossing of the resonance.
The excited beam dynamics becomes visible in particle tracking simulations, where the motion of particles shortly before their respective extraction is recorded (Fig. 5).As a consequence of the interplay between the excitation signal and the oscillatory particle motion-subject to amplitude and phase detuning-the periodic increase and decrease of the value H of the Kobayashi Hamiltonian is observed 32 .
The long coherence time of the sinusoidal excitation results in a much larger change of H on a short timescale compared to the diffusive nature of random kicks caused by an excitation signal with a bandwidth.As a result, the population of particles in the vicinity of the separatrix is greatly reduced, since particles coming too close to the separatrix will inevitably cross it, become unstable and get extracted.This depletion in the vicinity of the separatrix causes the extraction process to be less sensitive to ripples of the separatrix size, which greatly improves the spill quality.Table 2. Optimized signal parameters from the experimental investigation with tunes Q ex,i , bandwidths Q ex,i and period T FM .

Proposed excitation signal
Although the periodic in-and decrease of H in pure sinusoidal excitation is advantageous for the spill quality, it does not provide a sufficient net-increase for most particles to make for an efficient extraction.The extraction efficiency is defined as the number of of particles that are extracted during the simulation divided by the total number of particles simulated.In addition, the kick strength required to extract the beam using the sinusoidal excitation signal is an order of magnitude larger compared to the noise excitation.On the other hand, the diffusive nature of the noise excitation ensures a net energy transfer, providing a good extraction efficiency at an adequate signal level.To allow for an efficient extraction with reasonable kick strength while still improving the spill quality, a tradeoff has to be found.
We propose an excitation signal consisting of a band filtered uniform noise and two sinusoidal signals, in the following referred to as Noise++.For comparison of this new signal type, the parameters found in the experimental study (Table 2) are used in the simulation for the Noise++ signal.
Figure 6 shows that the Noise++ excitation achieves the same efficiency as the noise-only signal, i.e. it is able to the extract 95 % of the particles during the predetermined duration of the simulation.The missing 5 % are due to the feedback controller slowly increasing the signal level in the beginning of the simulation until the desired target rate is reached.Compared to other commonly used methods providing a net energy transfer, such as phase or frequency modulated signals, filtered noise has the advantage that there are no artificial periodic structures present in the signal which could introduce unwanted artificial low frequency ripples (compare "Experimental comparison" section).
The sinusoidal components reduce the particle density in the vicinity of the separatrix by the mechanism described above.While the effect is less strong in combination with the noise component, a significant reduction can still be achieved (Fig. 5).This improves the spill quality towards a Poisson-limited spill as shown in Fig. 6.The reduction of fluctuations in the spill rate is not limited to the time resolution of 50 µ s used in Figs. 2 and 6; it also manifests itself on smaller time scales as it reduces the probability of particles being extracted shortly after each other (pile-up).This aspect is highlighted in Fig. 7, where it is shown that the fraction of particles being extracted within sub-µ s time intervals is lower for Noise++ and sinusoidal excitation compared to the pure noise excitation.At multiples of f −1 ex a step-like increase is observed as a consequence of the excitation frequencies being imprinted onto the spill.For the considered average rate of 4 × 10 5 particles/s ≤ f ex this imprint does not affect the pile-up, such that for the case of sinusoidal excitation an almost Poisson-like distribution is achieved.
It was found in simulations and experiment, that the use of two sinusoidal components with slightly different fractional tunes enhances the effect compared to only a single one.To prevent beating, it is important to place the components in different betatron sidebands 23 .In practice, finding the exact frequencies that provide the best spill quality is an optimization problem.Therefore, adding more than two sinusoids makes the optimization lengthy without increasing the spill quality significantly.

Experimental comparison
The experimental comparison of excitation signals is carried out in July 2023 at the Cooler Synchrotron (COSY) in Jülich 33 .Table 1 lists the beam and optics parameters used for extraction of a coasting proton beam.A softwaredefined radio (SDR) based system and the GNU Radio software is used to generate the excitation signals 34,35 .This system also digitizes and records the signal of a plastic scintillator based particle detector in the extraction beamline.A built-in feedback (bandwidth approx.5 Hz) adjusts the excitation signal level in order to maintain an average spill rate of 5 × 10 5 p/s , which is sufficiently low for the used detector.
For each of the excitation methods presented in the following, the bound optimization by quadratic approximation (BOBYQA) algorithm 36,37 is used as a robust, derivative-free global optimizer in order to find the optimal excitation signal parameters.As objective function to be minimized, the coefficient of variation c v is used with t count = 500 µs and t evaluate = 4.2 s .The optimized parameters are listed in Table 2.
The output signal of the SDR is split and amplified with 42 dB 150 W power amplifiers.Inside a stripline unit, the resulting effective angular deflection of a traversing particle with a corresponding rigidity of Bρ = 3.3 Tm is 38 www.nature.com/scientificreports/with P SDR < 12.8 dBm being the signal power at the SDR output.

Sinusoidal excitation
As discussed in "Spill smoothing by coherent excitation" section, a purely sinusoidal excitation has a poor extraction efficiency and requires high signal levels.In the experiment, at least three sines are required to achieve extraction.Figure 8 shows the frequency spectrum of the resulting spill.Power supply ripples are prominently visible at 150, 300 and 600 Hz, originating from the three-phase bridge rectifiers commonly used for DC conversion 39 .Additional ripples appear at 2 and 3.6 kHz, which can be attributed to certain control loops of the COSY power converters.
The sinusoidal excitation strongly suppresses random fluctuations above about 500 Hz and also reduces the power supply ripples compared to noise excitation, confirming the spill smoothing as motivated by the simulation studies.However, it also introduces two artificial ripples at 11 and 26 kHz which are caused by the high signal levels required to extract particles, leading to clipping at the SDR's digital-to-analog converter (DAC) limit and saturation of the used RF amplifier.

Frequency modulated excitation
Frequency modulation (FM) is commonly used to overcome the efficiency limit of sinusoidal excitation by linearly sweeping the frequency in a sawtooth or zigzag pattern with period T FM 20,40 .However, the resulting periodic crossing of the machine tune leads to a pulsed extraction with strong artificial ripples at the sweeping frequency and its harmonics.
Dual FM aims to reduce this effect by using two carrier signals whose modulation is shifted by 180 • in phase 20 .Yet, the experiment shows that the sweeping frequency and harmonics are still imprinted onto the spill and dominate the spectrum (Fig. 9), even though the signal parameters were optimized to maximize the spill quality.While the addition of mono frequency sinusoidal components lowers statistical spill fluctuations and power supply ripples, it cannot counteract the artificial ripples introduced by the FM methods itself.

Phase modulated excitation
Phase modulation (PM) offers an alternative way of broadening the bandwidth of a sinusoidal signal.In its simplest form-random binary phase-shift keying (RBPSK) 41 -the phase is flipped randomly between 0 • and 180 • at fixed time intervals T flip = (�Q ex f rev ) −1 .Phase-shift keying is a standard encoding technique in com- munication system and therefore easy to implement.However, the phase flipping causes abrupt changes of the  extraction rate and introduces artificial ripples 23 .Figure 9 shows these ripples in the measured spill spectrum, coinciding with the first two harmonics of the flipping frequency for the optimized -3 dB bandwidth Q ex .
A variant with three carrier frequencies and reduced bandwidth has shown to improve the spill quality 22 .Considering the arguments in "Spill smoothing by coherent excitation" section, this improvement can be understood by the increased coherence time, representing a tradeoff between Poisson-limited sinusoidal excitation and an adequate extraction efficiency.
Figure 9 shows the improvement achieved using such a signal.The optimization algorithm converged to a parameter set where two carriers show a narrow bandwidth and one is broadband (Table 2).The reduced bandwidth can, however, not counteract the artificial ripples it causes in the first place, and even brings them to lower frequencies, where they are potentially more disturbing.

Noise excitation
SDR technology allows real-time generation of pseudo-random noise, providing an excitation signal free of periodic structures which could cause artificial ripples in the spill.A 3 rd order Butterworth infinite impulse response (IIR) filter limits the spectral power to the desired bandwidth.While free of artificial ripples, the incoherence of the noise causes a slow diffusion of particles towards the separatrix, as opposed to the fast approach when using coherent sinusoidal signals.This leads to a broad noise floor in the measured spill spectrum (Fig. 10) even after optimizing the signal parameters.
In analogy to the multi-RBPSK signals, the combination of three noise components at different sidebands helps to reduce these incoherent fluctuations and the influence of power supply ripples.
As shown in Fig. 10, the effect of sinusoidal spill smoothing is maximized with the Noise++ signal described in "Proposed excitation signal" section, as it suppresses external ripples while not introducing artificial lowfrequency ripples.

Comparison
Figure 11 shows the spill quality achieved for the signals discussed above.This allows to compare the achievable spill quality not only for the timescale at which the signal parameters were optimized ( t count = 500 µs ), but also globally independent of the choice of t count .Starting from the pure noise excitation as an incoherent signal, the spill quality improves with the RBPSK method, the 3 RBPSK with reduced bandwidth and the sinusoidal excitation (3 sine), as the signal becomes more coherent.The sinusoidal excitation, however, suffers from its high power requirements exceeding the hardware capabilities (Fig. 12) which leads to clipping and diminishes the spill quality as described in "Sinusoidal excitation" section.The dual FM method overcomes the power limitation and-having a long coherence time-is able to provide a better spill quality at t count = 500µs .However, the artificial ripples introduced by the sweeping result in a much worse quality on smaller timescales, where the sweeping period is resolved.The Noise++ method improves the spill quality on all timescales down to 5 µs by combining the sinusoidal spill smoothing with a method which does not introduce artificial low-frequency ripples.Its RMS power requirement is comparable to the classic RBPSK or dual FM methods (Fig. 12) and while the peak power required is slightly higher, it is still well within the margin of hardware capabilities.
While the presented excitation signals are carefully optimized, the accelerator optics and slow extraction process were not optimized excessively, leaving room for further improvement beyond the presented results, which are solely based on the tuning of the excitation signal.If sufficient excitation power is available, or the beam energy is sufficiently low, the simulation results suggest that pure sinusoidal excitation at multiple sidebands will provide the best possible spill quality, limited only by Poisson statistics.

Conclusion
The study of signals for RF-KO driven resonant slow extraction shows that sinusoidal signals lead to a modulation of the energy of betatron oscillations by means of the Kobayashi Hamiltonian.As a result, the phase space density in the vicinity of the separatrix is reduced, and the transition from stable to unbound motion by crossing the separatrix is sped up.This reduces the vulnerability of the extraction process to any external ripples modulating the separatrix size.It also avoids the low-frequency characteristic imprinted onto the spill as a consequence of the interaction of other excitation signals with the anharmonic motion of particles near the separatrix.Both effects improve the spill quality and reduce pile-up compared to excitation signals with a non-zero bandwidth.
The Noise++ excitation method is proposed to combine these benefits of sinusoidal excitation with a signal of non-zero bandwidth that ensures a good extraction efficiency, while not introducing additional artificial ripples at low frequencies.In this study, the first fair comparison of all commonly used excitation signals under identical machine conditions is performed with the aid of an automatic optimization algorithm.This experimental comparison demonstrates the capability of the new method to provide a smooth spill on all timescales.
The excitation frequencies Q ex ≈ n ± q x ( n ∈ N ) should generally be placed in different betatron sidebands with sufficient distance between them to prevent low frequency beating.Depending on the requirements, sidebands above a certain frequency should be chosen to prevent the time structure of the excitation-which in RF-KO generally is imprinted on the spill-from disturbing the downstream user of the beam.Once the sidebands are chosen, the signal parameters (fractional tunes and bandwidth) of the Noise++ excitation can be optimized for the respective accelerator optics.For this purpose, the SDR based tool with its built-in optimizer can be utilized as presented in this study 35 .
In addition to the simulations and experimental studies at COSY presented here, the authors have performed equivalent experiments at the Heidelberg Ion Beam Therapy Center (HIT)-including simulations-and at GSI SIS-18.These show qualitatively the same behaviour for the investigated sinusoidal spill smoothing and the Noise++ excitation signal, indicating that the presented findings are of general validity for RF-KO resonant slow extraction from synchrotrons.
Using a slightly higher excitation power than conventional methods, we expect the Noise++ excitation signal to provide the best overall performance in terms of spill fluctuation reduction in facilities using the RF-KO slow extraction method.When fluctuations on small timescales below 1 ms are not of concern, other signals like Dual FM achieve the same performance with less excitation power.The 3 RBPSK signal is particularly suitable for applications where excitation power is the limiting factor, and gives excellent results on timescales of 10 ms and above.

Figure 1 .
Figure 1.Schematic of resonant slow extraction.Left: Horizontal phase space at the electrostatic septum.Circles (filled for every 3rd turn) show the trace of a particle crossing the separatrix (red dashed), becoming unstable and jumping over the septum blade (black).Right: Aerial view of the beamline.The exciter controls the transition from the stable (thin line) to the unstable orbit (thick line), where the particle is deflected by the electrostatic (E) and magnetic (M) septa.

1 Figure 2 .
Figure 2. Spill rate (left) and spill spectrum (middle, right) from particle tracking simulation for comparison of noise and sinusoidal excitation.The indicated target rate is Ṅ = 5 × 10 5 /s .The zooms show the frequency region with the imprint of the respective excitation signal (see "Simulation framework" section).The dotted lines indicate the location of the resonance (red) and tune (black).

Figure 5 .
Figure 5. Normalized Hamiltonian distribution of particles before their respective extraction, simulated for the different excitation signals described in "Simulation framework" section.Noise++ is an intermediate excitation between noise and sinusoidal excitation (see "Proposed excitation signal" section).The average is plotted as white dashed line on top of each distribution, and in the right plot its time derivative in the close vicinity of the separatrix is shown.The separatrix is marked by the red dashed lines.

Figure 6 .
Figure 6.Comparison of spill quality for simulated noise, sinusoidal and Noise++ excitation.The quality metrics as defined in equation 3 are determined with t count = 50 µs and t evaluate = 50 ms .The number in brackets denotes the extraction efficiency, i.e. the percentage of particles extracted during the simulated spill duration.

Figure 7 .
Figure 7. Cumulative distribution function of delay between the extraction of consecutive particles.Results of simulated noise, Noise++ and sinusoidal excitations are shown.The dashed line indicates the ideal Poisson distribution.

Figure 8 .
Figure 8. Measured spill spectrum for sinusoidal excitation.Power supply ripples (filled triangles) and ripples due to saturation (outlined triangles) are indicated.

Figure 9 .
Figure 9. Measured spill spectra for frequency and phase modulated excitation.Power supply ripples (filled triangles) and ripples caused by the modulation (outlined triangles) are indicated.

Figure 10 .
Figure 10.Measured spill spectrum (left) and spill rate (right) for different noise-based excitation signals.Power supply ripples (filled triangles) are indicated.

Figure 11 .
Figure 11.Measured spill quality for different excitation methods at two distinct timescales (left) and as function of timescale (right).Signal parameters optimized for t count = 500 µs (vertical line).The respective Poisson limit are indicated.

Figure 12 .
Figure 12.Measured output power of the SDR for different excitation methods.Peak and RMS power are indicated alongside the distribution of signal levels throughout the spill (hatched violin plots).The corresponding deflection is given by equation 12.